Entanglement Entropy in the Two-Dimensional Random Transverse Field Ising Model 
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The scaling behavior of the entanglement entropy in the two-dimensional random transverse field 
Ising model is studied numerically through the strong disordered renormalization group method. 
We find that the leading term of the entanglement entropy always scales linearly with the block 
size. However, besides this area law contribution, we find a subleading logarithmic correction at 
the quantum critical point. This correction is discussed from the point of view of an underlying 
percolation transition, both at finite and at zero temperature. 
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The study of novel quantum phases and related quan- 
tum phase transitions (QPTs) is at the forefront of many 
recent developments in condensed matter physics. It re- 
lies heavily on the concept of entanglement entropy. 

A state |<J>) of a bipartite system A{J B is entangled if 
it cannot be described accurately in either subsystem A 
or B. A convenient measure of this entanglement is the 
entropy, S A = -Trp A \og 2 p A , where p A = Tr B |$)($|. 
Denote the linear dimensions of A [J B and A as M 
and L, respectively. An important question in quantum 
many-body systems is to study how S A (L) scales with 
L in the limit of M — > oo in different quantum phases. 
This question has been extensively investigated in one 
dimension (ID) [H, S [3, 0, i, @, 0] • There, it is now well 
understood that for non-critical systems, S(L) saturates 
to a constant as L — > oo; whereas in critical systems, a 
logarithmic modification stands out as the leading term: 
S(L) ~ InL, and its coefficient is associated with the cen- 
tral charge of the related (1 + 1) conformal field theory 
(CFT) [7]. In higher dimensions, it is generally believed 
that an area law holds at least for non-critical systems: 
the entanglement entropy scales as the area of the bound- 
ary between subsystem A and B, S(L) ~ This has 
been confirmed by studies on bosonic harmonic lattice 
systems [1, Q. For critical systems, the situation is more 
complicated. The area law is shown to be violated in free 
fermion systems with a finite Fermi surface [^, [Io| ■ But 
it still holds for fermionic systems without a finite Fermi 
surface [13], and critical bosonic systems [Til ] . 

In (i-dimensional system where the area law holds, 
S(L) ~ f s L d ~ x . Here f s is a boundary free energy deter- 
mined by the short-distance properties of the system, and 
is hence not universal. It is thus interesting to wonder 
about subleading [n| terms in S(L), where, maybe, uni- 
versal coefficients depending only on the model and the 
topological properties of the system could appear. For 
instance, it was recently found that in two-dimensional 
(2D) gapped systems, a subleading constant contribution 
in S(L) is related to the topological order [13(. Also, for a 
class of z — 2 conformal quantum critical systems in 2D, 
a universal logarithmic correction to the area law term 



has been found [3]. Clearly, the problem is not fully 
settled. 

It is of course also possible to investigate entanglement 
in quantum disordered systems. In a scries of studies 
in ID based on strong disorder renormalization group 



(SDRG) techniques [17j j . it was found that for the class 
of ID infinite randomness fixed points (IRFP), a InL 
term in S(L) is also present [3, EH. Similar results were 
also discovered for ID aperiodic systems flij ]. 

In this paper, we report on our study of the 2D random 
transverse field Ising (RTFI) model, and the numerical 
calculation of the entanglement entropy using the SDRG 
technique. The model is defined on a 2D square lattice 
with linear dimension M and open boundary condition. 
The subsystem A is a L x L square region located in the 
center of the square lattice. The Hamiltonian reads 



H 



(id) 



'hiSf. 



(1) 



The Ising coupling Jy and the transverse field hi take 
random values drawn from the following box shape dis- 
tributions: 



p(j) = e(j)-e(j-i), 

P(h) = ^[Q(h)-Q(h-h )]. 
ho 



(2) 



hase transi- 
Here the 



This model is known to have a quantum 
tion which is governed by an IRFP 1| 
critical point is tuned by ho- Starting from the origi- 
nal Hamiltonian Eq. [I] the SDRG finds the ground state 
by successively eliminating the highest energy degrees 
of freedom [H HJ. At each RG step, we look for the 
largest term in the Hamiltonian; its coupling (or field) 
is defined as the energy scale f2 at this step. If 57 = hi, 
the local spin is frozen in the eigenstate of Sf by the 
local field. It is then eliminated from the system, and 
an effective coupling Jj fc « max(Jjk, JjiJik/hi) is intro- 
duced between its two neighboring spins at sites j and 
k. If f2 = Jij, the two spins involved respond to the 
field uniformly, so that they are combined into a new 



2 



effective spin (or a ciuster). Then we effectiveiy elimi- 
nate one spin degree of freedom, and the local field at 
the new effective spin is h' { = hihj/Jij. These two dec- 
imation procedures are illustrated in Fig. [T] Numeri- 
cally the RG is processed until only one cluster is left in 
the system. The ground state then consists of indepen- 
dent clusters each of which is frozen into a GHZ state: 
\C(n)) = J= (| t)® n + | I)®") ■ Since each GHZ state 
will contribute either 1 (we take logarithms in base two) 
to the entropy if it consists of degrees of freedom in both 
subsystems A and B, or otherwise, calculating the en- 
tanglement entropy between two subsystems is reduced 
to a pure cluster counting problem: S(L) is proportional 
to the number of clusters N(L) that cross the bound- 
ary between the two subsystems. Finally S(L) is aver- 
aged over different disorder configurations. In practice, 
10 5 — 10 6 configurations are used. 
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FIG. 1: (color online) Basic RG transformations (see text for 
detail). Upper: energy scale is a field; lower: energy scale is 
an Ising coupling. 



The foregoing technique has been applied with success 
in the ID case. The technical difference in 2D is first 
that the quantum critical point is not exactly known. To 
locate it accurately we study the scaling behavior of the 
average magnetization m(M) [2l|. At the critical field 



m(2M) 



m(M) 



= 2~ 



(3) 



is independent of M, where x m is the anomalous dimen- 
sion of the bulk magnetization. Our result is given in 
Fig. [UJa). The critical field is estimated to be Hq = 
5.37 ± 0.03, with x m = 1.01 ± 0.05, which is consistent 
with a previous RG study 22]. The entropy S(L) is 
calculated for various values of hg. For both critical and 
noncritical ho, we find that the area law holds: S(L) ~ L 
in the leading term. The result of S(L)/L for different 
system sizes at critical /ig = 5.35 is shown in Fig.^b). 

This conclusion is quite different from the one in a 
recent study [23^ . where a double-logarithmic modifica- 
tion of the area law in the same model was reported at 
the critical point. We find that for small systems the 
double-logarithmic fit is reasonable, but that for system 



size M 128, S(L)/L increases definitely slower than 
lnlnL for L > 48, strongly suggesting S(L) ~ L in the 
limit of L — > oo, without modification, and that the ob- 
servation of Ref. [23[ is biased by finite size effects. We 
also note that our results are largely independent of the 
distribution of couplings, confirming the idea of universal 
behavior for S(L). 
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FIG. 2: (color online) (a): Finite-size scaling of magnetization 
ratio given in Eq. [3] (b): entropy per surface S(L)/L v.s. 
lnlnL at critical field = 5.35, The dashed line is a linear 
fit in In In L scale. 



Having established the validity of the area law in this 
system, it is natural to investigate subleading terms. We 
thus consider SS(L) = 2S(L) — S(2L), in which the terms 
linear in L cancel exactly. We find in both disordered and 
ordered phases that SS(L) saturates to a constant term, 
indicating S(L) — aL + c. Meanwhile, at the critical 
point, we find that SS(L) scales linearly as InL, suggest- 
ing 



S(L) = aL + blaL + c, 



(4) 



i.e., a logarithmic correction to the area law. The coef- 
ficient of this logarithmic correction is determined to be 
b = — 0.019±0.005 through finite-size scaling in Fig.[3Ja). 
To our knowledge, this is the first instance of such be- 
havior in disordered 2D systems, and the second instance 
in all 2D systems after the examples in Ref. 14] for a 
class of conformal quantum critical models with dynam- 
ical exponent z — 2. There is no reason to expect that 
the InL term we find in the critical RTFI model has 
much to do with the latter. This can be substantiated 
by calculating the amplitudes of the logarithmic term 
for different geometries, which obey some precise rela- 
tions in the case of Ref. |l4j. As an example, we con- 



sidered a cross shape geometry as shown in Fig. [3th). 
In this case as well, we can resolve a InL term in S(L) 
in addition to the area law contribution, with the co- 
efficient b cross = —0.08 ± 0.01. We can then calculate 
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and compare the ratios in our model, where we obtain 
across /bgquare ~ 4, and in the conformal quantum criti- 
cal models where b cross /b square = 3 exactly. This implies 
that the lnL term in S(L) in our model most probably 
has a different origin, a non surprising conclusion since, 
for the IRFP, z -> oo. 
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FIG. 3: (color online) (a): Scaling plots of 5S(L) to reveal 
the subleading term of S(L) at critical field /i§ = 5.35 for 
different geometries of subsystem A: a square (in (a)), and a 
cross shape (in (b)). Dashed lines are linear fits in InL scale 



with b s 



-0.019 and b c 



-0.08, respectively. 



To better understand the In L term in S(L) at the 
IRFP in 2D RTFI model, we notice that there is a strik- 
ing difference between the model in 2D and in ID. In 2D 
for any ho < hf, there is a finite-temperature phase tran- 
sition at T c (h ) [18]: the IRFP in 2D can then be consid- 
ered as an extension of this finite-temperature transition 
right down to T — 0. Through the SDRG, the transition 
to a ferromagnetically ordered phase can be mapped to a 
percolation transition in 2D [1 91 ] : the magnetic transition 
corresponds to the development of an infinite percolat- 
ing spin cluster during RG. It is widely expected that 
this percolation process at the IRFP (which occurs at 
energy scale fioo = 0) is different from the one at finite 
temperature, since at ho = h c the critical behavior is con- 
trolled by quantum fluctuations. This leads one to think 
of the IRFP as a type of "quantum percolation" , with 
fractal dimension df = 2 — x m « 1.0. For ho < /ig mean- 
while, the percolation takes place at finite energy scale 
Qoo ~ T c in the RG, and is expected to be in the univer- 
sality class of conventional classical percolation [r| . We 
have confirmed the classical percolation picture at finite 
energy scale by studying the scaling of largest active clus- 
ter size during RG. Some numerical results at ho = 3.2 
are presented in Fig. [4^a) . The percolation threshold is 
at Too = ln(f2o/^oo) = 5.27 ± 0.02, where the extracted 
exponents (3 and v take the values of classical percolation 
indeed. 



The number of clusters crossing the boundary between 
two subsystems can be investigated at these percolation 
transitions as well (even though it does not correspond 
to an entanglement entropy except when h — hg). It 
is easy to see then that a InL correction to the area 
law is expected and related to conformal invariance, even 
though it depends on more complicated parameters that 
the central charge and the topology (in contrast with the 
example in Ref. [H]). 
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FIG. 4: (color online) (a): Scaling of the largest active cluster 
size N c during RG shows a signature of classical percolation 
at a finite energy scale; (b): scaling of 8N at percolation 
threshold Too = 5.27 for ho = 3.2, A is a line interval on 
the boundary; (c): scaling of 8N at bond (upper) and site 
(lower) classical percolation thresholds, A takes the geometry 
of a square. In (b) and (c), dashed lines are linear fits in InL 
scale with corresponding slope b marked on the plot. 



To see this, consider the "baby" case where the sub- 
system A takes the geometry of a line interval of length L 
on the boundary of the lattice, so the boundary between 
the two subsystems is A itself. Now N(L) is simply equal 
to the number of clusters touching A. Its scaling can be 
studied using CFT techniques. For this, consider first the 
problem on the upper half complex plane, with A on the 
real axis. Use the well known expansion of the partition 
function of the critical Q-stete Potts model in terms of 
clusters, or equivalently, dense loops [24|: 



■Np 



(5) 



where every loop gets the same weight \[Q. Now 
introduce a boundary conditions changing operator 
(BCCO) [25| 4> y such that the two point function of tfi 
is defined through a sum similar to the one for Z, but 
loops touching the boundary between the two insertions 
Ci and (2 get a different weight y instead of ^/Q: 



(6) 



N P ((l,(2) 
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where N P ((i, £2) gives the number of loops touching the 
boundary between (1 and £2 located on the real axis. 
We expect the two point function to have the following 
scaling form: 

< MtiWvfo) >~e-'MICi-&l| Cl _^ J |-J%) i (7) 

where /(y) is the boundary free energy induced by the 
modified weight on the boundary, and the exponent h(y) 
is the anomalous dimension of the BCCO. Now differen- 
tiate the two point function of BCCOs with respect to 
the weight y, then take the limit y = \[Q. This leads to 



N b P (L) = aL + blnL, 



where 



-2y/Q 



dh(y) 



dy 



y=VQ 



(8) 



(9) 



and L = |Ci - Ca|- 

With the exact expression of h(y) [26| we obtain 



(10) 



For percolation (Q = 1), b perc = 4^ ~ 0.1378, repro- 
duces an early result by Cardy [27J . But in Eq. [TU] we 
generalize Cardy's result to general Q, and it is interest- 
ing to see that the In L term vanishes at Q = and Q = 4. 
It is also remarkable to see that b is related to the deriva- 
tive of the anomalous dimension of BCCO, but not the 
central charge. This In L term is also observed in our RG 
calculation at finite . In Fig. 0Jb) we show the scaling 
of SN = 2N(L) - N(2L) at r m = 5.27 for h = 3.2. b is 
estimated to be 0.15 ± 0.02, in agreement with the ana- 
lytical result. This further confirms that the universality 
class at finite is classical percolation. Interestingly, 
we find numerically for this case that \b\ < 0.01 at the 
quantum critical point, consistent again with the idea of 
a different universality class when — ► 00. 

Going back to the original problem where the subsys- 
tem A takes the geometry of a L x L square, we have 
not derived a similar analytical result for ordinary per- 
colation. But N(L) can of course be calculated numer- 
ically. To get better scaling for large systems, we turn 
to a direct study of percolation. In Fig. (He), SN data 
at percolation threshold are shown. As well expected, 
SN ~ InL is resolved for both bond and site percolation, 
and the coefficient of the In L term takes the same value 
b = —0.06 ± 0.01, in agreement with the idea that this 
term is universal. Note that we observe a negative value 
of b for subsystem A a square, just as in the case of the 
RTFI model. This is opposite to the sign of b in classical 
percolation when A is an interval. 

The observation of a In L term in N(L) in percolation 
makes the presence of a similar term at the IRFP most 
likely: there will in fact always be such a term at energy 



scale Ooq. When this scale is finite, the coefficient b takes 
the value of classical percolation, b = —0.06 ± 0.01. But 
when — ► 0, i.e., at the quantum critical point, quan- 
tum fluctuations become dominant, leading to a quan- 
tum percolation belonging to different universality class. 
A different b value, b = -0.019±0.005, reflecting this 
difference is then observed. 

In summary, we have calculated the entanglement en- 
tropy of a 2D RTFI model by using a numerical SDRG 
method. In contrast to what is claimed in a recent 
preprint, we find that the leading term of the entropy 
follows the area law and depends linearly on the block 
size L in both critical and non-critical phases. However, 
a In L correction to the area law is discovered at critical- 
ity While the presence of this correction may not have 
been expected from the entanglement point of view, it is 
very natural once the problem is reformulated geometri- 
cally. Indeed, the problem of counting clusters touching 
a boundary in 2D classical percolation is easily argued 
to give rise to sub logarithmic corrections, while the en- 
tanglement entropy in the RTFI model at criticality can 
be reformulated as a similar problem but in a different, 
"quantum percolation" universality class. 
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